BFKL pomeron in the external field of the nucleus in 

(2+l)-dimensional QCD 



M.A.Braun, A.N.Tarasov 
St. Petersburg State University, Russia 

^ ; March 3, 2013 

Abstract The behaviour of the pomeron propagator in the external nuclear field is studied 
in the (2+l)-dimensional QCD. It is shown that in the physically interesting case when the field 
does not vanish at large rapidities the propagator in the field vanishes much faster than in the 
vacuum, in agreement with the results found in the local Regge-Gribov model. However if the 
nuclear field vanishes at high rapidities the field does not change the behaviour of the pomeron 
^ ■ propagator. 

<D '■ 
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i — i; 1 Introduction 

In the framework of the perturbative QCD at small x and large N c the strong interaction is 
' realized by the exchange of hard BFKL pomerons which interact via the 3-pomeron vertex 

governing their splitting and fusion. One can sum all the corresponding fan diagrams in the 
quasi-classical approximation by means of the Balitski-Kovchegov (BK) equation for 7 (h)A 

■ scattering [1, 2] or by means of a pair of equations for AB scattering introduced in [3]. With 
the growth of energy the role of pomeron loops becomes important and one has to search for 
methods to take them into account. In our previous publications we noted that calculations of 
loops may become easier if one starts with the perturbative approach inside the nucleus from 

■ the start. At least in the local Regge-Gribov model the nuclear field transforms the supercritical 
^ pomeron with the intercept greater than unity into a subcritical one with the intercept smaller 

than unity [4]. As a result with the growth of energy the loop becomes relatively smaller, so 
that it can be calculated perturbatively. 

This result can be easily understood from the structure of the model in the nuclear sur- 
rounding. The Lagrangian of the Regge-Gribov model has a form 

£ = $t(JL + m$ + T/$t$($t + $) (1) 
oy 

where <3? and <&t are quantum fields which describe the pomeron, H is the free " Hamiltonian" 
and V - the 3-pomeron vertex. In the quasi-classical approximation one has <3?t = and <I> = $0, 
the latter given by the sum of all fan diagrams corresponding to the interaction with the nucleus. 
To calculate quantum correction one can shift the field — > <I> + $0 an d thus study the model 
in the nuclear background field ^o- This introduces a new interaction V$>o& corresponding to 
the annihilation of a pair of pomerons and, most important, the pomeron propagator must now 
be calculated in the background field <I>o- It turns out that this propagator, which in the vacuum 
exponentially grows with rapidity, in the background nuclear field exponentially vanishes. As 
a result in the background field pomeron loops do not grow with rapidity and can be easily 
controlled. 
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Of course the problem is to see if this result can be generalized to the perturbative QCD 
pomeron. Our numerical calculations have shown that the QCD pomeron propagator in the 
external field of the nucleus also vanishes at high rapidities in contrast to its behaviour in the 
vacuum [5]. Unfortunately the realistic QCD in 3+1 dimensions does not allow to confirm this 
result analytically, since both the relevant BK equation and the BFKL equation in the external 
field do not allow for analytic solutions. As was shown a some time ago the situation is improved 
in the (2+l)-dimensional QCD, where a solution to the BK equation can be obtained in a form 
which analytically shows its rapidity dependence [6]. 

In this note we use this known solution to study the behaviour of the BFKL propagator in 
the external field given by this solution and corresponding to the field inside a big nucleus. Our 
results show that this behaviour, as well as the behaviour of the nuclear field itself, very much 
depend on the initial conditions for evolution. If they are chosen in such a way that at large 
rapidities the fan-diagram-field tends to unity (the dipole S-matrix vanishes) then solutions of 
the BFKL equation in this field also vanish in this limit. However for different initial conditions, 
for which the fan diagram- field vanishes at large rapidities (the dipole S'-matrix tends to unity), 
this field does not change the behaviour of the solutions of the BFKL equation. A striking 
result of our study is that the interval of the coordinates in which the BFKL propagator in 
the external field is different from zero diminishes with the fall of energy. As a consequence 
pomerons cannot form loops in the nuclear field. This implies that the quasi-classical approach 
gives the complete solution to the quantum field theory of interacting pomerons in the nuclear 
field in 2+1 dimensions. 



2 Main equations 

As found in [6] the kernel of the BK equation (the rescaled triple pomeron vertex) is changed 
for the (2+l)-dimensional space according to 

n N r 2 

-i^ZJZT^ 2a.*c*(rSr - r 3 )6(r 3 - r™), (2) 

z/l '23' 31 

where r ™® x ( min ) [ s fa e max imal (minimal) of the coordinates r 2 and n; r 2 \ = r 2 — r \ etc. This 
vertex is a constant different from zero only when r 3 lies between r 2 and r±. This allows to write 
the BK equation in a simple form. Let S r2ri (y) be the S'-matrix for the interaction of a dipole 
stretched between spatial (1-dimensional) points r 2 > r\ at rapidity Y. The BK equation for S, 
which sums all the fan diagrams, is 

8S r 2 



dy 



J dr (s r2ro S rori - S r2ri ), y = 2iraY, (3) 

where Y is the rapidity. This equation can be further simplified if we introduce a function 

* r2ri (y) = e r2iy S r2ri , r 21 = r 2 - n > 0. (4) 

The equation for \I> is 



dr tf r2ro tf rori . (5) 



ri 



dy 

Considering ^ r2 r 1 (y) as a triangular matrix 

*ran(j/) =< r 2 \^{y)\ri >, r 2 > n, (6) 

we can rewrite Eq. (5) as 

* 2 (y), (7) 



dy 
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with a solution satisfying the initial condition 

^r 2 rM =Sr 2ri (y) = (0) (8) 

and given by the matrix formula 

*(y) = *(0)[l-y*(0)]- 1 . (9) 

To obtain the equation for the BFKL pomeron inside the nucleus we first rewrite Eq. (3) in 
terms of the scattering amplitude <& r2ri (y) defined by 

Snniv) = 1 -*ran(j/)- (10) 

We get 

u ^r 2 r\ I 



rr 2 

/ dr 



This scattering amplitude represents the sum of fan diagrams connecting the projectile with the 
nucleus target. This is precisely the field created by the nuclear background. 

Dropping the non- linear term one obtains the BFKL equation in 2+1 dimensions in the 
vacuum ^ 

= j dVQ^Pj-^TQ ~\~ Pron ~ P'r-2rij ■ (1^) 

The BFKL equation in the nuclear field is obtained when one adds to the right-hand side terms 
which correspond to the interaction of the pomeron with the nuclear field via the triple pomeron 
vertex _ _ r 

J dro(^P r , 2ro + .Fj-on -F > r . 2ri ^r 2 roPror\ — -Ff^ro^rori^ ■ (1^) 



BP f r l 



By Jri 

Expressing <E> via S we rewrite this equation as 
BP r r 2 



By 

Finally we introduce 



J dro^S r2ro P rori + Pr2r(,Sr () ri Przri^ • (1^) 
Qr2rM = e r2lV Pr2rMi (15) 



/ dr 



= {Q(y)Mv)}- (17) 



for which the equation simplifies to 

BQ r2 n _ 

By 

or in the matrix notation 

dQjy) 
By 

Solution to this equation with the initial condition 

Q(y,y') , = Q(y',y') 

y=y 

is found to be 

Q(y,y') = ^- 1 (y')^(y)Q(y',y')^(y)^- 1 (y > ) 

^T(y,y>)Q(y>y>)T(y,y>), (18) 

where we define triangular matrix 

T(y,y') = = *(y)*-Hy')- (19) 
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Note that ^(y) and commute at any y and y', since they both are functions of the same 

matrix *(0). 

In the following we shall consider a simple case when the nuclear 1-dimensional density is 
constant ("nuclear matter"). Then the initial BK function "I / r2ri (0) will depend only on the 
distance between the gluons 

\EV 2n (0) = VEV 21 (0). (20) 

As a result \EV 2n (y) will depend only on the distance p 2 \ at all values of y and consequently 
T r2ri (y,y') will also depend only on the distance r 2 \. As a matrix 

T r2ri (y,y') =6{r 21 )T r21 {y,y'). (21) 

3 The Green function of the BFKL equation (14) 

We shall first study the Green function G r2ri \ r > r i (y,y') which satisfies 

G r2ri | rW (y', y') = 5(r 2 - r' 2 )S(n - r[) (22) 

and we assume that r 2 = max{r2,ri} and r' 2 = maxjr^r^}. Note that this is not the pomeron 
propagator, which will be studied later. To find G we have to search for a solution of the 
equation in variables r 2 and r\ for P which satisfies (22). In terms of 

Fr 2ri (y,y') = e yV21 G r2ri y 2r/i (y,y') 

we have to find a solution of Eq. (14) which satisfies 

F r2ri (y',y') = e y' r ^9(r 21 )5(r 2 - r' 2 )5(n - r[). 

Obviously we should also have r' 2 > r[. 
According to (18) it is given by 

F r2ri (y,y') = ey'^e{r 21 )9(r' 2l ) J dr dr' T r2ro (y, y ')6(r - r' 2 )5{r' Q - r'^T^y^y') 

= ey'^9(r 21 )e(r' 21 )T r2r , (y, y ')T r[ri (y,y'). (23) 
The Green function itself will be given by 

G r2rilr > r[ (y,y') = e^i-^eCraije^jr^^y.y')^^^,^). (24) 

This expression is different from zero only in the interval 

Gr^y^iy^y') + only if n < r[ < r' 2 < r 2 . (25) 

Here r 2 and r' 2 are the maximal coordinates of the initial and final gluons and r[ and r' are 
their minimal coordinates. From (25) one concludes that the interval of coordinates where G is 
different from zero is restricted on both sides. This interval at y should be wholly inside the one 
at the initial rapidity y' . This remarkable property will have far-reaching consequences for the 
whole model, as will be discussed in the following. 

For the considered case of nuclear matter when T rari = T r21 the Green function G can be 
presented in the form which explicitly shows its symmetry in the initial and final gluons 

G r2riKr{ ( y , y ') = ey' s '-y s T S2 (y,y')T sl (y,y'). (26) 
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where 

s = max{r2,ri} — min{r2,ri} s' = max{r 2 ,r^} — min{r 2 ,r^}, 

s 2 = max{r 2 ,ri} - max{r 2 ,r' 1 }, si = min{r 2 , r^} - min{r 2 , ri} (27) 

and s, s', si and s 2 are all to be non-negative so that s 2 + si = s — s' > 0. 
At y' = we have 

T r2ri (y,0) = ([l-y^(0)]- 1 ) . (28) 

From the equation 

£ dr T r2ro (y,0)(6{r - n) - i/*ron(0)) = <5(r 2 - n) (29) 

we conclude that 

T r2ri (y,0) = 5(r 2 - n ) + Y r2ri (y), (30) 

where Y(y) is a smooth function of r 2 and r\. Presence of the 5-function allows to get the 
correct right-hand side, which otherwise would be zero at r 2 = r\. Using the fact that the BK 
solution depends only on the difference r 2 i we find 

T r21 (y) = 5{r 2l ) + Y r21 {y). (31) 

To find Y r {y) we set up an evolution equation in y. Obviously 

(Ul iM ' i}) = (*(0)T 2 (y,0)) r = jf" dn* r _ ri (0) jH dr 2 T ri _ r2 (y,0)T r2 (y,0), (32) 



tf r (0) + 2 [ dr^r-rMYrM + f dr^ r - r M f 1 dr 2 Y ri _ r2 {y)Y r2 (y) (33) 

JO JO JO 



dy 

or in terms of Y 

dY r (y) /•'■ /•' /■'' 

<9y 

with the initial condition 

y r (0) = 0. (34) 

From this evolution equation one can numerically find Y r (y) at any y by the standard Runge- 
Kutta procedure. At arbitrary y and y' 

T(y,y') = [l-y / ^(0)][l-y*(0)]- 1 . (35) 

From this we find 

T r (y,y , )=5(r)+X r (y,y'), (36) 
where X(y,y') is a smooth function linear in y'\ 

X r (y,y') = Y r (y)-y'h r (0) + f dr'y r - r ,Y r ,(y)\ . (37) 
L Jo J 

In absence of the interaction with the nucleus S'-matrix corresponding to the BK equation 

turns into unity. In terms of the triangular matrices it implies 

sW 1 (y) = 9(r 2 -r 1 ) (38) 

and does not depend on y. From this we find in the vacuum 

T^( y , y ') = 6(r) + (y-y')ey r e(r). (39) 

Using this and (24) we obtain the Green function in the vacuum: 

G^ rilr , r ,(y,y') = e-(y-y'y^{5(r 22 ,)5(r vl ) 

+ (y ~ y')5(r vl )9{r 22 ,) + (y - y')5{r 22 ,)9(r vi ) + (y - y') 2 e{r 22 ,)9{r vi )), (40) 

where it is assumed that r 22 ', ry\ > 0. This expression coincides with the one previously obtained 
in [6]. As expected in absence of the external field the Green function depends only on the 
difference y — y' . 
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4 Pomeron propagator in the nuclear field 

In contrast with the Green function G defined by the equation 



(—+h)g=1, (41) 



where H is the BFKL Hamiltonian in the external field, the pomeron propagator g is defined by 

(— + H)g = pi 2 P 2 2 , (42) 



so that it is related to the Green function considered previously as 

g(r 2 , n \r' 2 , r[) = V{ 2 V 2 2 G(r 2 , n \r' 2 ,r[) = G(r 2 , n \r' 2 , r[)V^ 2 V 2 2 , (43) 
or equivalently 

V?V!<7(r 2 ,n|r£,ri) = 5 (r 2 , n^, ^) V 2 2 = G(r 2 ,n|r 2 ,ri), (44) 

where in the last equality the derivatives act on the left. Obviously relations (44) determine g 
non- uniquely. One can add to it any function linear in r 2 — r' 2 or in r\ — r[. Our choice will 
consist in selecting g to be different from zero in the same intervals of initial and final gluon 
coordinates as G, namely s, s', s\, s 2 > where s are defined by (27). As we shall find the 
propagator g defined in this way at y > y' vanishes at large separation of points. Linear terms 
that can be added to it and different from zero outside the interval (25) violate this condition 
and seem to us unacceptable on formal grounds, since correspond to strong infrared singularities 
at all values of y. 

To make our choice we notice that in the one-dimensional transverse coordinate space func- 
tion h(r) = V~ 2 5(r) can be presented as 

h(r) = r9(r). (45) 

Indeed we have 

V 2 /i(r) = h"(r) = (rd(r))" = rS'{r) + 28{r) = 5{r). 

This implies that the Green function g(y, r 2 ,r\\y', r 2 ,r^) may be chosen to satisfy the BFKL 
equation 

(— + H)g{y, r 2 ,n\y', r 2 ,r[) = (r 2 - r' 2 )6{r 2 - r' 2 )(r[ - ri )e(r[ - n), (46) 

where it is assumed r 2 i, r' 21 > 0. Note the inverse order of the initial and final coordinates in 
the second factor on the right. It is this order that guarantees the desired property of g to be 
different from zero only in the interval (25). 

We can solve this equation in the same manner as when searching for the Green function G. 
We introduce matrix 

f r2 rMv') = 0(r 21 ) e y^g r2rilr;A (y,y') (47) 
with r 2 i, r' 21 > 0. It will satisfy Eq. (46) with the initial condition 

WiA V') = 0(r 21 )ey' r ^9(r' 21 )(r 2 - r' 2 )6(r 2 - r' 2 )(r[ - n)0(ri - n). (48) 

According to Eq. (18) the solution is given by 

fr 2ri (y,y') = E 0(r 00/ )T r2ro (y,yy o2/ 9(r 02 y io 9(r[ )ey^ (49) 

ro,r' n 
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where we use the notation r o' = tq — r' and r' 10 = r[ — r' . 

Taking into account that T r2ri depends only on the distance r 2 \ between the points we find 

gr^r'^y') = e(r 2l )9{r' 2l )e- r ^y-y"> jdrj dr'6(r 21 - r - r>)6(r 22 , - r)9{r vi - r') 

e-y'^ +r '\r 22l -r)(r vl -r')T r {y,y')T rl {y,y'). (50) 
Integration over r and r' is restricted by the last two ^-functions to the interval 

< r < r 22 >, O <r' < r vi . 

From this we conclude 

r + r'< r 22 , + r V \ = r 21 - r' 21 < r 21 
so that 9(r 2 \ — r — r') = 1. As a result this expression factorizes 

9r 2riKr[ (y,y') = e- r ^y-y^U r . 22 ,(y,y')U rill (y,y'), (51) 
where we have defined a matrix 

U r {y,y') = ^ dr'(r-r')T r ,(y,y>)e-y' r '. (52) 
Jo 

In the notation symmetric in the initial and final gluons this can be rewritten as 

9r^\ri 2 r^y,y') = e- s ^U sl (y,y')U S2 (y,y'). (53) 

where s, s', s\ and s 2 are defined by Eqs. (27). 
We have seen that 

T(y,y') = [l-y'^(0)][l-y^(0)}- 1 (54) 
and that the second factor contains a (^-contribution 

[1-y^iO)}' 1 = l + Y(y). (55) 

From this we find 

T(y,y') = l + X(y,y'), X(y, y') = Y(y) - y'tf(0) - y'*(0)Y(y), (56) 
so that T r2Tl contains a 5 term. If T r , 2Tl depends only on the distance r 2 \ we find 

T r (y,y')=5(r)+X r (y,y'), (57) 
where X r {y,y') is a smooth function of r linear in y'\ 

X r (y,y') = Y r (y) - y'h r (0) + [ dr>^ r „ r ,Y r ,(y)\. (58) 
L Jo J 

Separating the 5 function from T one gets 

U r (y,y')=r+ f dr' X r ,{y,y')(r - r')e^ V . (59) 
Jo 

In absence of the nuclear field, using (39) one obtains 

e {y-y')r _ i 

u?\y-y') = e y _ yl - ■ (60) 

From this using (51) we find the vacuum propagator, which naturally depends only on the 
difference y — y', which we denote as simply y 



g r2riKr[ (y) = l -{e-y^+e-y^-e-y^-e-y^'}, (61) 
where it is assumed that r\ < r[ < r' 2 < r 2 . 



1e-25 1 ' ' 1 1 1 

2 4 6 8 10 

r 

Figure 1: S-matrix S r (y) evolved by the BK equation from the initial function Eq. (64) with 
7 = 0.1 as a function of r. Curves from top to bottom correspond to y = 1, 3, 5, 7 and 9 

4.1 Numerical results 

Both the nuclear field characterized by the BK S'-matrix S r2ri (y) and the pomeron propagator 
are determined by the coupling to the nucleus contained in the BK S'-matrix at y = 0. In 
the nuclear matter on physical grounds at r = one requires S r (0) r= o = 1 corresponding to 
vanishing of integration for zero-dimensional dipoles. As noted in [6] the behaviour of S with the 
growth of energy crucially depends on the behaviour of S r (0) S r (0) r _>o ~ 1 — cr 7 . In particular 
for 7 < 1 it was found that the asymptotical behavior at y » 1 was 

S r (y) ~ e-ciXi+T^ 1 -^. ( 62 ) 

On the other hand for 7 > 1 the S-matrix was expected to tend to unity in the large region 
r << 3/T- 1 . 

In our calculations we have taken S r (0) in the form 

S r (0) = e~ rl (63) 

Numerical results show that for 7 < 1 the asymptotical behaviour (62) is well fulfilled. In Fig. 

1 we show the behavour S r (y) with 7 = 0.1 for values of rescaled rapidity y = 1, 3, 5, 7 and 9. 
In fact they are well described by (62) practically at all values of r and y. However for 7 > 1 
we find the asymptotical behaviour somewhat different from what was expected in [6]. In Fig. 

2 we present S r (y) with 7 = 2 at the same rapidities y = 1, 3, 5, 7 and 9. They all fall onto the 
exponential dependence on r 

S r (r)~e~^ r . (64) 

In contrast to the case 7 < 1 the slope b(y) falls with y but not to zero but to the finite value 
~ 0.33. The same picture is observed with still higher values of 7. 

Passing to the pomeron propagator we note that it depends on three independent variables. 
Assuming that 77 < r[ < r' 2 < r 2 , they can be taken as any three of the four differences r 2 i, r' 2 i, 
r 2 2' and ryi, taking non-negative values and constrained by the relation 



r'21 + r 2 2' +rn = r 2i . 



(65) 




To see both r— and y— dependence we shall first study integrals of g at y' = over r 2 2' and 
r±'i with fixed r 2 i = r 

gi (y,r) = e~y r f d Tl U ri (y, 0) f * dr 2 U r2 (y, 0) (66) 

JO JO 

and with fixed r 2 i = r 

rrm-r rr rn -r-r 1 

g 2 (y,r)) = e-y r d ri U ri (y,0)e~^ dr 2 U r2 (y, 0) e -^ 2 , (67) 

JO JO 

where r m is the upper limit of the integration, (formally infinite, but finite in the course of 
calculations and taken r m = 20). 

Our numerical results for g^ and g^ with 7 = 0.1 and y = 1, 3, 5, 7 and 9 are shown in Figs. 
3 and 4 respectively. Similar results with 7 = 2.0 are illustrated in Figs. 5 and 6. To compare 
in Fig. 7 and 8 we plot the same quantities for the vacuum case. In all cases the integrated 
propagator diminishes with energy and coordinate. 

To clearly see the influence of the field and the role of value of 7 in Fig. 9 we show values 
of g(y,0) integrated over all three independent variables from (65). To facilitate the comparison 
results with 7 = 0.1 and 7 = 2.0 were rescaled to coincide with the vacuum result at y = 1. We 
observe a clear distinction between the cases 7 < 1 and 7 > 1 previously stressed in [6]. With 
7 < 1 we find that the propagator in the nuclear field indeed goes to zero much faster than in 
the vacuum as advocated in our previous studies in the 3+1 dimensional case. In contrast with 
7 > 1 we do not observe any influence of the field. This goes in line with the assertion in [6] 
that this case corresponds to asymptotic vanishing of the dipole interaction. 

5 Loops 

In 3+1 dimensions in the theory of interacting BFKL pomerons the interaction is given by 
Sj = ^[ dy d2nd2r 2 2d 2 r " rl2 NM^(y, ri ,r 2 )} ■ fty, r 2 , r 3 )<p(y, r 3 , n) + (h.c). (68) 

7T J r 23 r 31 
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1e-25 u ' ' 1 1 1 

2 4 6 8 10 

r 

Figure 3: Pomeron propagator integrated over the final gluon coordinates gi(y,r) as a function 
of r. The initial BK function Eq. (64) is taken with 7 = 0.1. Curves from top to bottom 
correspond to y = 1, 3, 5, 7 and 9 
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Figure 4: Pomeron propagator integrated over the initial gluon coordinates #2(2/, r) as a function 
of r. The initial BK function Eq. (64) is taken with 7 = 0.1. Curves from top to bottom 
correspond to y = 1, 3, 5, 7 and 9 
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Figure 5: Same as Fig. 3 with 7 = 2.0 in the initial BK function Eq. (64) Curves from top to 
bottom correspond to y = 1, 3, 5, 7 and 9 
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Figure 6: Same as Fig. 4 with 7 = 2.0 in the initial BK function Eq. (64) Curves from top to 
bottom correspond to y = 1, 3, 5, 7 and 9 
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Figure 7: Same as Fig. 3 for the vacuum case. Curves from top to bottom correspond to y 
3, 5, 7 and 9 




Figure 8: Same as Fig. 4 for the vacuum case. Curves from top to bottom correspond to y 
3, 5, 7 and 9 
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Figure 9: Fully integrated propagator g(y, 0). The two curves with 7 = 0.1 (the lower one) and 
7 = 2.0 (one of the two practically coinciding upper curves) are rescaled to have the same value 
with the vacuum curve (the other of the two upper curves) at y = 1 



Accordingly the pomeron self-mass S is given by 
^{y,r2,n\y ,r 2 , ri ) = -[ J / —5—5 g —g{y,r 2 , r 3 \y , r 2 , r 3 )g{y, n, r 3 |y , r l5 r 3 ) 

If J r 23 r 31 r 2'3' r 3'l' 

(69) 

(the minus sign reflects the absorptive character of the triple pomeron interaction) . 

Passing to 2+1 dimensions we substitute the vertex function rf 2 /^23 r i3 by the difference of 
9 functions, as found in [6] and indicated in (2). The pomeron self-mass becomes 

/•max{r2,ri} rmax{r' 2 ,r ^ } 

Z(y,r 2 ,r 1 \y',r 2 ,r' 1 ) = -(8wN c a 2 s ) 2 / dr 3 / dr 3 

Jmin{r2,ri} J mm-fr^r^} 

5(y, r 2 , r 3 |y', r' 2 ,r' 3 )g(y, r 3 , n r 3 , ri) (70) 

It is trivial to see that this expression is zero due to the properties of the pomeron propagator. 
In fact take r 2 > T\ and r 2 > r[. The first propagator is different from zero only at r 3 > r 3 . and 
the second propagator is different from zero only at r 3 > r 3 . So their product is zero. 

This circumstance is obvious since the interval where the propagator is different from zero 
is extending on both sides with the growth of y as schematically shown in Fig. 10, a. Splitting 
the pomeron in two generates the picture shown in Fig. 10.6. It is evident that after splitting 
as the rapidity diminishes the two new pomerons can never interact with each other (points r 3 
and r 3 ' in Fig. 10,6 cannot coincide). So loops do not exist in our theory. 

6 Conclusions 

In 2+1 dimensions solutions of both the BK equation describing the nuclear field and BFKL 
equation in this field can be obtained in the closed form, expressed via infinite-dimensional 
matrices on the points dividing the coordinate interval. This makes it possible to find both the 
Green function of the BFKL equation and the pomeron propagator. It is important that the 
initial condition for the propagator at y = y' in principle admits addition of terms linear in 
one of the two final coordinates. We have proposed a choice which in our opinion is physically 
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Figure 10: Schematic regions in the (y, r) plane where the pomeron propagator (a) or a pair of 
propagators after splitting (b) are different from zero 

motivated, since it has no infrared singularities at y — y' > 0. With this choice calculation of the 
propagator shows that the influence of the nuclear field very much depends on the behaviour of 
this field. If at large rapidities the nuclear S'-matrix tends to zero (black disk limit) the field 
indeed makes the propagator fall with rapidity much faster than in the vacuum. On the other 
hand if the nuclear S'-matrix tends to unity at high energies (zero interaction with the nucleus) 
the field in fact does not change the behaviour of the propagator in comparison to the nuclear 
case. 

An important consequence of our choice of the propagator is that the region where it is 
different from zero does not change with energy and restricted to the interval of final coordinates 
which contains inside the one of the initial coordinates. As a result pomerons cannot form loops 
of any sort. Thus with this choice the quasi-classical solution of the BK equation gives the 
complete solution of the dipole interaction with the nucleus. The question if other choices of 
the propagator are admissible from some points of view remains open and will be considered in 
future studies. 
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